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ABSTRACT 

Measurement of the time delay between multiple images of a gravitational lens 
system is potentially an accurate method of determining the Hubble constant over 
cosmological distances. One of the most promising candidates for an application of 
this technique is the system B0218+357 which was found in the Jodrell Bank/VLA 
Astrometric Survey (JVAS). This system consists of two images of a compact radio 
source, separated by 335 milliarcsec, and an Einstein ring which can provide a strong 
constraint on the mass distribution in the lens. We present here the results of a three- 
month VLA monitoring campaign at two frequencies. The data are of high quality and 
both images show clear variations in total flux density, percentage polarization and 
polarization position angle at both frequencies. The time delay between the variations 
in the two images has been calculated using a chi-squared minimization to be 10.5±0.4 
days at 95 per cent confidence, with the error being derived from Monte-Carlo simu¬ 
lations of the light curves. Although mass modelling of the system is at a preliminary 
stage, taking the lensing galaxy to be a singular isothermal ellipsoid and using the new 
value for the time delay gives a value for the Hubble constant of 69l|^g kms“^ Mpc“^, 
again at 95 per cent confidence. 

Key words: gravitational lensing - cosmology: observations - cosmology: miscella¬ 
neous - quasars: individual: B0218-I-357 


1 INTRODUCTION 

Long before the first gravitat ional lens was discovered, it 
had been shown (Refsdal 1964) that measurement of a time 


delay between the images of a lens could be used to calcu¬ 
late the Hubble constant, Hq, independently of any other 
distance determination to the lens or lensed object. This 
technique has so far been applied predominantly to the Dou¬ 
ble Quasar B0957-I-561 (Walsh, Carswell & Weymann 1979) 
where a long-running controversy as to t he length of the time 
delay has only recen tly been resolved ( Kundic et al. 1997 
Haarsma et al. 199S). In the first of these two papers, an Hq 
of 64±13 kms“^ Mpc“^ was derived based on an optically 
determined time delay of 417±3 days (at 95 per cent con¬ 
fidence). However, the deflecting mass is complicated (com¬ 
prising a galaxy and a galaxy cluster) and modelling it sat¬ 
isfactorily has proved difficult and constitutes the biggest 
source of error on the value of Hq derived from this system 
at this time. 

Systems containing an Einstein ring are ideal candidates 
for determining Ho as the presence of the ring can firmly c on- 
strain the mass distribution in the lens (Kochanek 199C). A 


good example of this is B0218-I-357 (Patnaik et al. 1993), 
first identified as a gravitational lens through observations 
carried out as p art of the Jodrell B ank/VLA Astrometric 
Survey (JVAS) (Patnaik et al. 1992). This lens system has 
a simple morphology (see Fig. which consists of two com¬ 
pact images (A and B) of a strongly variable flat-spectrum 
radio core and a steep-spectrum Einstein ring, the diame¬ 
ter of which is the same as the separation of the compact 
components, 335 milliarcsec (mas). This is the smallest sepa¬ 
ration yet found in a galactic-mass gravitational lens system 
and as a consequence, the time delay between the variations 
in components A and B is small. The ring is believed to 
be an image of part of the extended structure of the kpc- 
scale radio jet and will therefore vary on much longer time 
scales than the variations seen in the images of the compact 
cores. The deflecting mass comprises a single isolated galaxy 
which, in contrast to B0957-I-561, can be modelled relatively 
simply. The galaxy is also almost certainly a spiral because 
radio absorption observations have shown that the column 
density of absorbing material is very high (Carilli, Rupen & 
Yanny 1993, Wiklind & Combes 1995) and because a large 
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Figure 1. VLA 15 GHz radio map of B0218+357. As well as the 
two compact components (A to the right) and the Einstein ring, 
also clearly visible is a (non-lensed) radio jet to the south. 


differential Faraday rota tion measure exists between A and 
B at radio wavelengths (Patnaik et al. 199c). The redshifts 


of the le ased object and leasing gala xy are well determ ined 
at 0.96 (Lawrence 199C) and 0.6847 (Browne et al. 1993) re¬ 
spectively. As these are relatively low for a lens system, the 
assumed cosmology introduces less uncertainty in the value 

of Hq than with oth er systems. _ 

Previous work (Corbett et al. 1996) derived a time de¬ 


lay of : 2 days (B lagging A) with a la error of ±3 days 
from VLA observations ot the percentage polarized flux at 
15 GHz. In this paper we present new results of a three- 
month monitoring campaign conducted with the VLA at 
8.4 and 15 GHz. 


2 OBSERVATIONS AND DATA REDUCTION 

B0218+357 was observed with the VLA in A configuration 
between the months of October 1996 and January 1997. 
Observations were taken at two frequencies, 15 GHz and 
8.4 GHz, each with a bandwidth of 50 MHz. Each band is 
further split into two intermediate frequencies, or IFs, sep¬ 
arated slightly in frequency. With resolutions of 120 mas at 
15 GHz and 200 mas at 8.4 GHz, the variable components 
A and B are easily resolved and can be monitored for varia¬ 
tions in total flux density, percentage linear polarization and 
polarization position angle. In all, data were obtained at 47 
epochs, with an average spacing between observations of ~2 
days. 

The observing strategy was kept as consistent as pos¬ 
sible over the period of the monitoring. Each epoch con¬ 
tained two observations of 3C84, one before and one after 
B0218+357, for amplitude and phase calibration purposes. 
The 3G84 data are also used to correct for instrumental 
polarization on the assumption that 3C84 is unpolarized. 
3C119 was observed as a ‘control source’ to check for varia¬ 
tions in the flux density of 3C84 (each epoch is calibrated as¬ 
suming the same flux for 3C84) and to calibrate the polariza¬ 
tion position angle of B0218+357. 3C119 is a steep-spectrum 
source known to contain a very weak core (Hen-dong et al 


be very small. An innovation compared to the previous mon¬ 
itoring campaign was that antenna pointing offset observa¬ 
tions were made of each source to ensure that gain variations 
associated with pointing errors were minimized. B0218+357 
was observed for ~20 minutes at 15 GHz and ~3 minutes at 
8.4 GHz resulting in expected rms noise levels of 0.19 mjy 
and 0.12 mJy respectively for one IF. 

Calibration was performed using the NRAO Astronom¬ 
ical Image Processing Software package AlPS and each IF 
was calibrated separately. As the primary amplitude and 
phase calibrator, 3C84, is slightly resolved at the observed 
frequencies, the first step in the calibration process was to 
make the best possible map of 3C84. This was done by per¬ 
forming a first order calibration of several epochs of 3C84 
data and combining these to produce a self-calibrated map 
of the source. The clean components so derived were then 
used as a model input to the AiPS task calib to derive tele¬ 
scope gain solutions at each epoch of observation. These gain 
solutions were then applied to each of the target sources, 
3C119 and B0218+357. The calibration of the 15 GHz total 
flux density data was improved significantly by correcting 
for the gain of the VLA antennas changing with elevation. 
At lower frequencies the gain/elevation dependence is much 
weaker and at 8.4 GHz a correction made little or no differ¬ 
ence to the gain solutions. No correction at this frequency 
has therefore been applied. 

The flux densities of components A and B of B0218+357 
were calculated by fitting a model of the source direc tly to 
the calibrated uv visibility data using the difmap (Shep¬ 


herd 1997) model fitter (for Stokes I) and its AIPS equiva¬ 


lent UVFIT (Stokes Q and U). The model was kept as simple 
as possible and consisted of two point sources (separated 
by 335 mas at a position angle of 67 degrees) plus a broad 
Gaussian to represent the ring. The model fitting consisted 
of varying only the flux densities of these three components. 
Whereas the positions of components A and B are known 
to very high accuracy from VLBI observations, the size and 
position of the ring at each frequency was determined by 
performing an initial fit to each epoch’s data and allowing 
the ring parameters to float. As no systematic variations in 
the ring parameters were seen, the size and position of the 
broad Gaussian were then fixed at the average value derived 
from all the epochs. 

Visually the model provides a good fit to the data at all 
epochs, especially at 15 GHz as the ring flux density is low 
and the ring heavily resolved on most of the baselines. The 
per degree of freedom (x^) of fits at this frequency were 
typically less than 2. At 8.4 GHz the ring emission is stronger 
relative to components A and B and less separated from the 
point sources, resulting in poorer fits (x^ ~ 3-7). In the case 
of 3C119 the source can be reasonably approximated by a 
single Gaussian. 

At this point we average the results found for the two 
IFs. Bad epochs could be identified by a particularly high x^ 
or through inconsistent ring or control source flux densities 
and/or position angles. In most cases, the problem was then 
traced back to the original data and rectified. Because of ei¬ 
ther bad weather or instrumental problems, eleven 15 GHz 
and eight 8.4 GHz epochs of total flux density measurements 
have been removed from the time delay analysis. A partic¬ 
ular problem occurred with the pointing offset observations 
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for the first three epochs which has rendered the total flux 
density data unusable. However, the polarization informa¬ 
tion was obtainable for all but one of the epochs lost in 
total flux by amplitude self-calibrating the data, the results 
of which seem compatible with the other epochs. This gives 
45 epochs of polarization data at 15 GHz and 44 at 8.4 GHz. 


3 THE RADIO ‘LIGHT CURVES’ 

The total flux density, percentage polarization and polar¬ 
ization position a^le light curves for components A and B 
are shown in Fig. (15 GHz) and Fig. ^ (8.4 GHz). As well 
as the light curves themselves we discuss in some detail the 
way in which the error bars on the points in each light curve 
are derived since these in the end determine the confidence 
we will have in the derived time delay. 


3.1 Total flux density 

The variability and time delay signature are well illustrated 
by the 15 GHz total flux density light curves. Fig. a) and 
b). The main trend in total flux density (seen at both fre¬ 
quencies) is a rise of about 10 per cent followed by a short 
plateau stage, and a sharp decline which is only just caught 
in the B component data. Also, prominent in the 15 GHz 
A component light curve are two dips in the flux density 
around days 20 and 30, both of which show up clearly in the 
B component data, but delayed by about 11 days relative 
to A. These features can be seen to a lesser extent in the 
8.4 GHz light curves. Fig. ^ a) and b). 

The variations in total flux density seen in the light 
curves for B0218-I-357 should be compared with those for 
the control source 3G119 shown in Fig. ^ a) and b). For 
both frequencies the rms scatter is .^1 per cent compared 
to the ^10 per cent changes seen in B0218-I-357. Fig. ^ 
appears to show that the total flux density of 3C119 in¬ 
creases by ~1 per cent over the period of the observa¬ 
tions. We do not believe this is a real change in 3C119, 
but an artifact arising from a gradual decrease in the flux 
density of 3C84 which was used as the flux density cali¬ 
bration source. The total flux density of 3C84 has in fact 
been falling almost continuously since about 1983 (Univer- 
sity of Michigan Interactive Radio Observatory Database , 
http://www.astro.lsa.umich.edu/obs/radiotel/umrao.html) 
and this manifests itself as the apparent small rise in the 
flux density of 3C119 over the period of the monitoring ob¬ 
servations. We tried normalising the data by dividing the 
B0218-I-357 total flux densities at each epoch by the corre¬ 
sponding 3C119 total flux density, but this had the effect of 
slightly increasing the scatter in the B0218-I-357 total flux 
densities and so has not been applied. 

The fact that normalizing the B0218-I-357 data by the 
3C119 light curve does not produce an improvement sug¬ 
gests that the residual scatter in the 3C119 data is due to ob¬ 
servational, calibration and model fitting uncertainties and 
is therefore a good measure of the flux density measurement 
errors, both for 3C119 and B0218-1-357. We therefore use the 
scatter on the 3G119 light curve to calculate an error bar for 
the total flux density measurements at 15 GHz. If we fit a 
straight line to the 3C119 data, the scatter around this is 18 
mjy or 0.9 per cent. As calibration errors are multiplicative. 


the error bar for each epoch in Fig. g a) and b) has been 
set equal to 0.9 per cent of the flux at that epoch, typical 
values of which are 8 mJy for component A and 2 mJy for 
component B. 

Using the same method to obtain the error bars at 
8.4 GHz does not work well. This is because the model fit 
to the B0218-I-357 data is so much poorer (as measured by 
yf) than that to 3C119. This indicates that model fitting 
errors dominate over calibration errors and the scatter in 
the 3G119 data at this frequency (11 mJy or 0.4 per cent) 
is far too small to be a good estimate of the errors on the 
B0218-I-357 8.4 GHz flux density points. In the end, we be¬ 
lieve the chi-squared minimization technique itself that we 
use to calculate the delay between components A and B (de¬ 
scribed in Section 4.1) gives the best estimate of the errors. 
This is because we expect a of approximately unity if 
the errors on the data points are correctly chosen. For a 
‘moderately’ good fit, the statistic has a mean of v and 
a standard deviati on of y/Tu where v is the number of de¬ 
grees of freedom (e.g.. Press et al. 1992). Therefore to be 
conservative we choose our errors such that at b^st fit 
is one standard deviation below v. For the B0218-I-357 to¬ 
tal flux density data set which contains 38 epochs and is 
fitted with two parameters (delay and flux density ratio), 
this gives a required x^ equal to 0.76 which corresponds to 
errors of approximately 7 mJy for component A and 2 mJy 
for component B. 


3.2 Percentage polarization 

The 15 GHz polarization curves also contain variations over 
short time scales which provide an independent constraint 
on the time delay. The time series of percentage polarization 
at 15 GHz broadly consists of two features that stand out 
from the noise. Seen in component B, Fig. d), the first is 
a rise and fall of greater than 1 per cent, taking place over 
a period of about 20 days, shortly after the beginning of the 
monitoring observations. The second is of similar magni¬ 
tude and duration to the first and can be seen most clearly 
in Gomponent A between days 70 and 90, Fig. ^ c). This 
corresponds to the sharp drop in total flux density seen at 
both frequencies. The 8.4 GHz percentage polarization re¬ 
sults, Fig. ^ c) and d), are of a much poorer quality than 
those seen at 15 GHz, due to the combination of a poorer 
Gaussian model fit to the data and lower polarized flux den¬ 
sity. Although there are some small variations that can be 
discerned, they are of such low significance compared to the 
other data sets that the 8.4 GHz polarization light curves 
have not been used in the main analysis below. 

The percentage polarization light curves of 3C119 at 
15 GHz and 8.4 GHz are shown in Fig. ^ c) and d). Clearly 
the polarization calibration procedure works well though 
there are some residual systematic trends of ~0.1 per cent. 
This suggests that perhaps the assumption that 3C84 has 
zero and non-variable polarization breaks down and that 
our derived polarization residuals are being contaminated 
by very small intra-day variations in 3C84. However, this 
will have a negligible effect on the B0218-I-357 polarization 
data as the variations in this source are much larger and oc¬ 
cur at different times and on different time scales than those 
seen in the 3C119 light curves. 

We use the rms difference between the measurements of 
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Figure 2. 15 GHz light curves. Top row: Total flux density of B0218+357, a) component A, b) component B. Middle row: Percentage 
polarization, c) component A, d) comnone nt B. Bottom row: Polarization position angle, e.) component A. f^ co mponent B. (The data 
plotted in Figs ^ and Q can be found at [ittp: //multivac.jb.man.ac.uk:8000/ceres/data_from_papers/0218.html ) 
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Figure 3. 8.4 GHz light curves. Top row: Total flux density of B0218+357, a) component A, b) component B. Middle row: Percentage 
polarization, c) component A, d) component B. Bottom row: Polarization position angle, e) component A, f) component B. 
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Figure 4. Top row: Total flux density of the control source 3C119, a) 15 GHz, b) 8.4 GHz. Bottom row: Percentage polarization of 
3C119, c) 15 GHz, d) 8.4 GHz. Note that we do not show the polarization position angle curves since 3C119 itself was used as the 
calibrator for the position angle of B0218+357. 


percentage polarizations obtained in the independent IFs in 
order to estimate the error bars for each data point. These 
differences will arise from thermal noise and model fitting 
errors and will, therefore, represent a lower limit to the error 
bar for each component. We also correct each IF difference 
for a constant offset that exists between the two IFs, the 
origin of which is unclear. At 15 GHz this offset is equal to 
0.02 per cent for component A and 0.07 per cent for com¬ 
ponent B. At this same frequency the estimated errors are 
equal to 0.05 per cent for component A and 0.17 per cent 
for component B. 


3.3 Polarization position angle 

By far the most prominent feature of the variation in po¬ 
larization position angle at 15 GHz is the gradual rotation 


between days 35 and 80 seen in both component A and B, 
Fig. ^ e) and f), covering a range of about 10°. As with per¬ 
centage polarization at 8.4 GHz, the polarization position 
angle variations at 8.4 GHz, Fig. H e) and f), are poorly de¬ 
fined compared to those at 15 GHz and are therefore also 
excluded from the time delay analysis. As the polarization 
position angle of B0218-I-357 is calibrated by that of 3C119 
we do not show plots of the 3G119 polarization position an¬ 
gle. 

The error bars on the 15 GHz and 8.4 GHz polarization 
position angle data are calculated in the same manner as 
that described in Section 3.2 for the percentage polarization 
data. Error bars on the 15 GHz data are 0.58° and 0.80° for 
components A and B respectively. 
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Figure 5. against delay for total flux density, 15 GHz. 


4 TIME DELAY ANALYSIS 


4.1 The chi-squared method 

In order to quantify the time delay between components A 
and B, a chi-squared analysis has been performed on each 
data set. Each point in one light curve, A{ti), is shifted by 
some delay r (between —10 and 40 days in steps of 0.1 of a 
day) and paired with a linearly interpolated value from the 
other, interpolation being required due to the uneven 

and sparse sampling of data. All points for which ti + r > fjv 
and ti — T < ti, where N is the total number of epochs, are 
excluded. As there is no reason to prefer interpolation in 
one time series over the other, the analysis is then repeated, 
but by shifting B{ti) by —r, interpolating in A{ti) and then 
averaging the results of the two passes. Since the measured 
chi-squared is inversely proportional to the number of over¬ 
lapping data points, delays for which the overlap is small 
are weighted down relative to those for which the overlap is 
large. Furthermore, as well as shifting the data temporally, 
the data are also scaled along the y-axis in order to produce 
the best alignment of the A and B light curves. Scalings are 
determined by hrst calculating the ratio (total flux density 
and percentage polarization) or difference (polarization po¬ 
sition angle) between each A and B pair in the unshifted 
data set and averaging over the entire light curve. A range 
of scalings that bracket this hgure are then determined and 
sampled in steps of 0.01 for total flux density and percent¬ 
age polarization and steps of 0.1 for polarization position 
angle. For each trial value of the delay, the value of this 
scaling is found which minimizes chi-squared, so calculat¬ 
ing the flux density ratio/de-polarization/differential rota¬ 
tion between A and B. An example plot of against delay 
is shown in Fig. ^ for total flux density, 15 GHz. 

As well as the chi-squared analysis described above we 
have also used a variety of other methods to determine the 
time delay, all of which have been used in the past for the 
analysis of B0957-I-561 monitoring data. We have calculated 
the cross-correlation function of the A and B light curve 
data , the discrete correlation function (Edelson fc Krolik 


1986) and the dispersion measure (Pelt et al. 1996), the 


last two specihcally taking into account the fact that the 
A and B time series consist of unevenly sampled data and 


hence do not require interpolation. These techniques all give 
results consistent with those found using a chi-squared min¬ 
imization and so in this paper we choose to concentrate on 
the results found from the chi-squared analysis. 


4.2 Error on the delay - Monte-Carlo simulations 

It is important to have an objective assessment of the conh- 
dence limits on the derived time delays. To do this we have 
adopted a Monte-Carlo approach taking into account both 
the uncertainty introduced by the measurement errors and 
those introduced by the gaps in the sampling of the light 
curves. We do this independently for each data set since the 
errors are different in each case. 

The first step is to try to create the best estimate of the 
true light curve of the object by combining the A and B light 
curves, making use of our estimate of the time delay and by 
scaling with our estimate of the image flux density ratio, de¬ 
polarization fraction or position angle difference. Examples 
of combined light curves are shown in Fig. 1^. The second 
stage is to re-sample with replacement^ the light curves at 
intervals chosen randomly from the distribution of real sam¬ 
pling intervals used for the original VLA observations. We 
choose the new sampling intervals so that the total length of 
the new light curve does not exceed that of the original. We 
also constrain its total length to be within 2.5 days of the 
original. Thus we ensure that all simulated light curves are 
of comparable length. Using these sampling intervals, the 
combined light curve is re-sampled twice, once with no de¬ 
lay and once using the delay estimated from the real data. 
Note that, by re-sampling the combined light curve with¬ 
out any smoothing, i.e. using linear interpolation, we take 
into account automatically much of the spread in the real 
data introduced by the effects of calibration errors etc. The 
resulting light curve is then de-magnihed, re-polarized or re¬ 
rotated by the appropriate amount used originally to create 
the combined light curve. New realizations of the A and B 
light curves are thus created. The samples forming the light 
curves for each epoch are then perturbed by adding random 
noise drawn from a normal distribution with a certain stan¬ 
dard deviation. This standard deviation is equal to the error 
bars on the points, the derivations of which are described in 
Section 3. From this perturbed set of light curves the best-ht 
delay and scaling is then calculated in exactly the same way 
as for the real data. This process is repeated many times, 
every time re-sampling the simulated data set. In this way 
frequency histograms are produced which should correspond 
to the error distribution of the measured delay. If anything 
this approach over-estimates the uncertainty in the time de¬ 
lay since there is some ‘double counting’ of the errors; they 
contribute to the spread in the combined light curve which is 
then resampled, and again when the resampled light curves 
are perturbed according to the errors on the points. 

Frequency histograms for the delay data are shown in 
Fig. Each is the result of 5,000 realizations and is plotted 
with a bin width of 0.2 day. The distributions for both de¬ 
lay and magnitude scaling are highly non-Gaussian and so 
cannot be described completely by their full width at half 


* A sampling interval once chosen is returned to the distribution. 
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Figure 6. Frequency histograms of delay (5,000 realizations) for a) Total flux density, 15 GHz, b) Percentage polarization, 15 GHz, c) 
Polarization position angle, 15 GHz and d) Total flux density, 8.4 GHz. Delays are binned in 0.2 day intervals. 


Table 1. Best-fit delays and magnitude scalings with associated 68.3 per cent confidence limits from the chi-squared analysis described 
in Section 4.1 and the Monte-Carlo simulations described in Section 4.2, for all data sets. ^A/B, ^A—B. 
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Figure 7. against delay for the simultaneous chi-squared min¬ 
imization. 


maximum, which would correspond to a Icr confidence inter¬ 
val. Note that each histogram has been plotted on the same 
axes thus clearly illustrating the relative merits of each data 
set. 

The time delays, magnitude scalings, associated 68.3 
per cent confidence limits and for each data set are shown 
in Table where the delay or magnitude scaling is equal to 
the median of the Monte-Carlo distribution and the con¬ 
fidence limits are calculated by counting inwards from the 
tails of the distributions until each tail contains 15.85 per 
cent of the total distribution. 

These results show that in general our derived time de¬ 
lays are compatible at la. Although the percentage polariza¬ 
tion and position angle data sets are not formally compatible 
at this confidence level, both are individually consistent with 
the time delays determined from the other data sets. 


4.3 The final delay estimate 

In order to arrive at a final best estimate of the time de¬ 
lay, the results from each of the four light curves (total flux 
density, 15 GHz and 8.4 GHz, percentage polarization and 
position angle, 15 GHz) are combined. We utilize a ‘simulta¬ 
neous chi-squared minimization’ technique which makes use 
of all the information contained in the against delay plots 
by summing the x^ found at each delay for each dataset. 
The minimum of the resulting distribution (see Fig. is 
that which is most consistent with all data sets and which 
represents our best-fit delay. This is equal to 10.5 days. 

To establish the error on this delay, we perform a 
Monte-Carlo simulation similar to that described in Sec¬ 
tion 4.2, forming a combined light curve from each of the four 
light curves, but using the best-fit delay found in the simul¬ 
taneous minimization for each. These are then re-sampled 
as before to produce realizations of the original light curves 
and another delay found by simultaneously minimizing chi- 
squared. This is repeated for a total of 5,000 realizations. 
The result from this simulation is plotted as a frequency 
histogram in Fig. ^ and the confidence limits given in Ta¬ 
ble At 95 per cent confidence the error on the time delay 
between images A and B in B0218-I-357 is ±0.4 days. 


Table 2. Allowed ranges in the combined delay for different con¬ 
fidence levels - simultaneous method. 


Confidence (per cent) 

Delay (days) 

68.3 

10.3 - 10.7 

90 

10.2 - 10.8 

95 

10.1 - 10.9 

99 

10.0 - 11.1 


o 



Figure 8. Frequency histogram illustrating the results of the 
simultaneous chi-squared minimization. Delays are binned in 0.1 
day intervals. 


As a cross-check on the simultaneous chi-squared 
method we take a weighted average of the 5,000 delays which 
comprise the results of the individual Monte-Carlo simula¬ 
tions (plotted in Fig. The delays are averaged realization 
by realization, weighting each by the inverse square of the 
full width of the Icr error bar (see Table Though the de¬ 
lay distributions obtained using this method are skewed to 
higher delays, the medians of the distributions using both 
methods are compatible at Icr. Our conclusion is therefore 
that at 95 per cent confidence the best estimate of the delay 
is 10.5 ± 0.4 days. 


5 DISCUSSION AND CONCLUSIONS 


Using new VLA observations of B0218±357, a value for the 
time delay between the variations of the two components 
in the gravitational lens system has been measured that is 
consistent with, and which represents an order of magni- 
t ude improvement on , the previous estimate of 12 ± 3 days 
( [Corbett et al. 1991: ). However, there is still room for im¬ 
provement and observations covering a longer time period 
and with more frequent sampling than those presented in 
this work will allow the error to be further reduced. 

Using the best estimate of the delay with the best-fit 
magnitude scalings, composite light curves are created by 
delaying each component A epoch by 10.5 days (see Fig. 

It is clear that this delay produces a good match between 
the A and B light curves, especially in the case of total flux 
density and polarization position angle at 15 GHz. An exam- 
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Figure 9. Combined light curves, a) Total flux density, 15 GHz, b) Percentage polarization, 15 GHz, c) Polarization position angle, 
15 GHz, d) Total flux density, 8.4 GHz. As in Figs ^ and g, component A and B measurements are represented by open and fllled squares 
respectively. 


ination of the percentage polarization composite light curve 
shows the match here to be poorer and suggests that the 


As further constraints we use the 15 GHz VLBI observa¬ 


tions of the ma s structure of the images (Patnaik, Porcas 


error b irs in this case have been under-estimated. This is _ fc Browne 1995). Both components A and B consist of two 


supported by the relatively high at best fit of the per¬ 
centage polarization light curves (see Table and could 
explain the formal Icr discrepancy between the time delays 
derived from this data set and from the polarization position 
angle data set. 

To obtain an estimate of the Hubble constant, we model 
B0218-I-357 using the singular isothermal ellipsoid (SIE) 
mass model, as described in Kormann, Schneider & Bartel- 
mann (1994). As constraints on the mass model we use the 
8.4 GHz flux density ratio A/B of 3.57 (Table |^, assuming 
an error of 0.015, which is approximately the 68.3 per cent 
confidence range. (Using the ratio 3.73 obtained from the 
15 GHz data increases the predicted delay by ~2 per cent.) 


subcomponents (A1,A2 and B1,B2) and all four are resolved. 
The flux density ratio is assumed to be the same for Al/Bl 
and A2/B2. We also make use of the positions of the sub¬ 
components as well as the deconvolved sizes and position 
angles of Al and Bl. Relative positional accuracies for the 
subcomponents of 0.1 mas are assumed. The centre of mass 
of the lens is a free parameter in the model. 

We create 5000 ‘data sets’, adding Gaussian-distributed 
errors to the observed positions, subcomponent sizes and 
flux density ratios. Using these artificial data sets, we solve 
for the minimum solution of the mass model parameters 
and the time delay between Al and Bl. The delay error 
distribution is shown in Fig. ^ and the resulting median 
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Figure 10. Delay distribution tunction obtained tor the lens 
model piU'timeLei'a Haled lii Tiible - 


ring. Given this caveat, however, combining the model re¬ 
sults with the observed delay of 10.5 ±0.4 days gives a value 
for the Hubble constant of 69ilijg kms“^ Mpc“^, where the 
errors are the formal 95 per cent confidence limits. This is 
for Ho = 1, Ao = 0 and r; = 1 (completely homogeneous uni¬ 
verse). Since the source and lens redshifts are relatively low, 
even generous departures from these fiducial values have a 
negligible effect on the above value for Hq. 

The above value is consistent with those found from 


other lens systems, e.g. B0957±561 (Kundic et al. 1997) 
B1608±656 (Fassnacht et al., Koopmans fc Fassnacht, 


in preparation) and PKS 1830-211 (Lovell et al. 1996) 


all of which lie within the range of approximately 60- 
70kms“^ Mpc“^. The situation with the remaining lens sys¬ 
tem for which Hq has been measured, PG 1115±080, is un¬ 
clear as several models have been published wi th a large 
spread (42-80) in the corresponding va lues of Hq (Schechtei 


et al. 1997; Keeton & Kochanek 1997). The value from the 


B0218±357 system is clearly capable of considerable im¬ 
provement, mainly in refinement of the lens mass model, but 
also by increasing the accuracy of the delay determination. 


Table 3. Median values of mass model parameters with associ¬ 
ated 68.3 per cent confidence limits. Lens centre coordinates are 
given relative to component A. Position angle is measured east 
from north. 


Parameter Value 


Lens centre {x) 
Lens centre (y) 
Velocity dispersion 
Axis ratio 
Position angle 


-9 

=;+4 


115_g mas 

no.rtlf kms-l 


r» y7-j-0.08 
'- 0.11 
qo+9 
-19 


-48° 


values of the mass model parameters are listed in Table ^ 
The errors indicate the range containing 68.3 per cent of the 
parameter distribution function. 

The centre of mass obtained from the model co¬ 
incides neither with the centre of the radio ring nor 
with the galaxy centre derived from an HST NIG- 
MOS observation (GASTLeS gra vitational lens database, 
http://pluto.harvard.edu/castles/). We do not regard these 
‘offsets’ very seriously at present. An offset of the model cen¬ 
tre from the centroid of the radio ring emission could arise 
if, in the source plane, the extended emission giving rise to 
the ring is asymmetrically placed with respect to the dia- 
mond caustic. An example of this is B1938±666 (King et 
al. 199^. Since existing radio maps do not define the ring 
structure very clearly it is difficult to tell if the ring is sym¬ 
metric from the radio images. The galaxy centre obtained 
from the HST NICMOS image is also hard to interpret in 
the lig ht of possible blending of the galaxy cor e with the B 
image (Xanthopoulos, Jackson & Browne 1998). 

The model predicts a median delay of 7.2'ti2'Q/h days 
between Al and Bl; Hq is equal to lOOh kms“^ Mpc“^. The 
quoted errors indicate the range containing 95 per cent of 
the delay distribution function as shown in Fig. ^ They are 
formal errors in the sense that they are based on the assump¬ 
tion that the effect of the lens is correctly modelled by the 
singular isothermal ellipsoid potential, an assumption which 
will be testable in the future when further modelling con¬ 
straints are obtained from better radio maps of the Einstein 
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